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Abstract 


In this paper we present a new and efficient method by combining pseudo 
differential operators and Haar wavelet to solve the linear and nonlinear dif- 
ferential equations. The present method performs extremely well in terms of 
efficiency and simplicity. 
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1 Introduction 


In this paper we have offered a method based on Haar wavelet in combination 
with pseudo differential operators to solve ODEs. 

Recently, wavelet methods have proved useful as widely applied tools for 
the computation of numerical approximations. Different types of wavelets 
and approximating functions have been proposed for the numerical solution 
of differential equations using various methods [2,8),9,[5120], such as Galerkin 
and collocations methods, the filter bank methods [20], and the Kalman fil- 
ters. In [§], the Haar wavelet method was applied for solving the boundary 
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and initial value differential equations. Lepik in [14,4] applied Haar wavelet 
to solve differential equations. Also, he used Haar wevelet for solving evolu- 
tion and integral equations in [1,[L3] and PDEs in [1]. 


pseudo differential operators are considered as natural extensions of lin- 
ear partial differential operators. The study of pseudo differential operators 
grows out of research in 1960s on the singular integral operators. The pseudo 
differential operator is the generalization of linear partial differential opera- 
tor, with roots entwined deep in solving differential equations in the classic 
and quantum mechanics [4,5]. 


At first by applying Fourier inversion formula to the given differential 
equation we obtain its associated pseudo differential equation. Then in the 
obtained Pseudo differential operator we expand f (f is assumed to be the 
solution of ODE) into Haar series, and we employ the collocation method 
for the equation to get a numerical solution for ODE. By some numerical 
examples the advantages of our tool in comparison with Haar wavelet method 
are shown. 


The paper is organized in 5 sections in the following way. In section 2 
we recall some preliminaries. Some numerical solutions for some examples 
of linear and nonlinear differential equations by using our new method are 
presented in sections 3 and 4, respectively. Finally error analysis is given in 
section 5. 


2 Preliminaries 


2.1 Wavelets 


Haar wavelets are the simplest wavelets among various types of wavelets. The 
Haar wavelet family on [0, 1] is defined as follows: 


1 for x € [&1, 2), 


h(x) = —1 for DE [€2, 3], 
0 elsewhere, 
where : 
k k+5 k+1 
&=—, = 2 =, 
m m m 


and m = 2/, 7 = 0,1,...,J, indicate the level of the wavelet, and k = 
0,1,...,m-—1 is the translation parameter. The index 7 is calculated accord- 
ing to the formula i = m+k+1. The maximal value of i is 2M where M = 2/. 
It is assumed that the value i = 1 corresponds to the scaling function of Haar 
wavelet h1(a) = 1 for all x € [0, 1]. 
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2.2 Fourier Transform and Pseudo Differential 


Operators 


One of the leading ideas in the theory of pseudo differential operators is to 
reduce the study of properties of a linear differential operator 


P=) ca(z)0, (1) 


which is a polynomial in the derivatives 0, with constants c, depending on 
x, to its symbol 


p(a,€) = > cale)(i€)*, 


which is a polynomial in the phase variable € € R with constants depending 
on the space variable «. 


Fourier transform is the main tool in the theory of pseudo differential op- 
erators. Let f be integrable function; the Fourier transform and the inverse 
Fourier transform of f are defined as follows: 


fQ) =3U© = I eE fan) 
and 
Sflla) = 5 fe fleas, (2) 


where z,€ € R. According to the inversion formula we have 


FGA) =F] =F (3) 


Let P be the linear differential operator in (M) using the inversion formula, 


PI = Yeale)Oty. fe Fleds 


Seal) 55 ie (ig) "el § fl@)ae 


= 2 [ c&p(a,e) (de. 


2n Jr 


Here p is considered as a symbol, and 


(Pla, Da)A(o) = fe (x, fede 
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defines the associated pseudo differential operator. For more details one can 
refer to [I]. 


2.3 Combination of Haar Wavelet and Pseudo 
Differential Operators 


Let f € L?[0, 1); then by (2) and (8), it can be approximated by 
fla) = 5 fe FOUe, (4) 


since f € L? then f € L?, and for level J we can write 


2M 


f= S- Aphy(E) 


where a, = i hy(€)f(€)dé, p=1,...,2M. By replacing this value in (@) we 
obtain 


2M 1 
= iw 
fa)= Yo [et holeas. (5) 
With differentiation of (G) we have 
2M 2M 1 
eo [ ibe hp(6)dé, f(z) =) _ ap / (i€)?e'*5 hy (€)d€, 
p=1 . 


in other words 
1 * 
= > a i; (i£)"e" hy (E)dE, forall neN. (6) 
p=1 o 


Now consider the linear differential equation Sv” _, Aa(x) f(a), by replac- 


ing the obtained value for f, f’,..., f in this equation we obtain 
n 2M 
So Aa(o) f(a Saf (x, Ee" hp (dé, 
a=0 
where p(a,€) = 33" _) Aa(x)(i€)*. This new form of pseudo differential op- 


erators is our main tool in order to obtain numerical solution for linear and 
nonlinear differential equations in the next sections. 
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2.4 Solving Linear Differential Equations 


Consider the differential equation of order n 
Yo Am (2)y (x) = F(), (7) 
m=0 


where x € [0,1] and y("-Y(0) = c,,...,y(0) = cp. The interval [0,1] is 
divided into 2M subintervals with the length Ax = sy 
According to the equality (@) we have 


2M 1 
yl (2) = Soap | (g)"etMag 6), (8) 
p=1 ” 
where ap, with p = 1,...,2M are unknown. Performing integration on (8) 
from 0 to a, yields 
2M 1 ; 1 
yea) = Lal | (yeh (@)ag— f (i€)"*Ap(E)dE] + cn (9) 
p=1 


continuing the integration of (G), we get 


2M 1 1 
uf (a) = Yo agl f(igyeMh(g)ag— ar" Fie) Ayla 
p=l1 
—m-—1 1 
D etyer™ (0) (10) 
for m = 0,1,...,n — 1. Putting the values obtained for y™, y"-),...,y 
and the collocation points 
(1 — 0.5) 
= ~—____ 1=1,2,...,2M 11 
TZ] 2M ) ae) ’ 5) ( ) 


in equation (M1), the system below is obtained. 


2M n 
Soa » Sin (2X) = f (x1) 
p=1 m=0 


for! =1,...,2M, where 
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Smn() = Am (2) | (ig) hy(€)de — 0-1 | (i€)"" hp (€)a€]. 


By calculating the coefficients ap, p = 1,..., 214, and replacing them in (10), 
for m = 0, an approximation numerical solution is obtained. The matrix 
form of the system (12) is 


(32 Wn) =F. 
m=0 
where W,,, signifies 2M x 2M matrices as follows: 
Wn (p,l) = Am(21)Sm(p, !), =O, Lyon 
such that 
Sm(p,!) = Sm(21), 
Sacpst) = f GeyeShgl eas, 


and F and a are 2M—dimensional vectors 


n—-1ln-v-1 


FW) =f + LO GletvOA(eo 
v=0 a=0 : 


a= (a1,...,@2m). 


Also, the matrix representation of (I), for m = 0, is 


Y= aSo + T, 
where the vector T is 
n-1 Tl 
T() => —(a)%y¥(0), 1 =1,..., 2M. 
rere. 


Remark. For calculating the integrals, for example, 
1 . 
[Wee hy dg, b= 0,15... 
0 


we can use of the points 7 = c», (=1,...,2M@ +1. In other words, 


2M 


1 Ts+1 : 
| (ig)FetEn,(€) = SO / (i€)Pe™Ehy (€)dE. 
0 s=1"7s 
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Since for € € (Ts, 7s41) we have h,(€) = hp(xs); then 


1 2M Hie . 
| (i)*e"* hp (€) = D7 hp (2s) / (ig)Fet™ de, 
" s=1 Ts 


3 The Numerical Solution of Linear Differential 
Equations 


In this section, some examples of linear differential equations are solved by 
using the new method. To get the error estimates, problems with the known 
exact solution ye,(x) are considered. 


Example 1. Consider the differential equation 
(x — 1)y!(x) + (x? + 3)y(x) = x4 + 7x? — 22 + 6, (13) 
where y(0) = 0. This equation has the exact solution yez(x) = 24+ 2?. 


Now, we seek solutions for the equation (1) based on our hybrid method. 
In this equation 


Wo(p, l) = (x7 + 3)So(p, l), 
Wi(p,l) = (a1 — 1)Si(p, 0), 
F(l) = af + 7x? — 2a, + 6. 


The obtained results and absolute errors, for some values of J = 2, are shown 
in Table 1 and Figure 1, respectively. 


Example 2. Consider the differential equation 


(x — 1)y""(x) + (w@ + 1)y'(®) + y(x) = (22? + 40 — 1)e®, 
(14) 


This equation has the exact solution yex(#) = xe. In this equation 
Wo(p,!) = So(p,!), 
Wi (p,1) = (a1 + 1)S1(p,!), 
W2(p, 1) = (a1 — 1)S2(p,!), 


F(l) = (2x? + 4a, — 1)e”' — (1+ 227). 
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Figure 1: The absolute errors for example 1. 
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Figure 2: The absolute errors for example 2. 
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Table 1: Haar solution, Exact solution, The new method solution for equation (2) 
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(a/16) Haar solution Exact solution The new method solution 
1 2.0087906524 2.0039062500 2.0039062497 
3 2.04280445967 2.0351562500 2.0351562500 
5) 2.11078811462 2.0976562500 2.0976562499 
7 2.21748292022 2.1914062500 2.1914062500 
9 2.38353209120 2.3164062500 2.3164062499 
11 2.79651329536 2.4726562500 2.4726562500 
13 0.16626536933 2.6601562500 2.6601562496 
15 4.80205593409 2.8789062500 2.8789062511 


The obtained results and absolute errors, for some values of J, are shown in 
Table 2 and Figure 2, respectively. 


4 Nonlinear Differential Equations 


In this section we adopt the new approach for nonlinear differential equations. 
A numerical solution of nonlinear differential equations can be obtained by 
using the new hybrid method. We consider a nonlinear form of differential 
equations as 

K(z,u(2),u'(2),... + u™ (2) = f(x) (15) 
such that a € [0,1] and u("-)(0) = ep,...,u(0) = cg. K is nonlinear respect 
to one of the functions u, u’,... uv and the function f is given. Replacing 
x by 2; in the equation (1), for 1 = 1,2,...,2M, gives 


K (a1, u(2,), u'(a1),... + u™ (a21)) = f(x), 


where 21, ¥,..., and x,y are the collocation points introduced in (MH). 
Now, by applying (8) and (@) for the functions u, u’,...,u', we get a system 
with 2M equations for calculating the coefficients ap. In general, this system 
is nonlinear, and we use the following approach to solve it. 
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Table 2: Haar solution, Exact solution, The new method solution for equation () 


(a/16) Haar solution Exact solution The new method solution 
1 0.0668034 0.0665309 0.0665309 
3 0.2271089 0.2261681 0.2261682 
5) 0.4290282 0.4271368 0.4271369 
if 0.6809506 0.6776132 0.6776133 
9 0.9929719 0.9872182 0.9872183 
11 1.3777507 1.3672570 1.3672570 
13 1.8556966 1.8309970 1.8309970 
15 2.4355610 2.3939901 2.3939901 


We assume that (14) is previously solved for the level J — 1, where the 
number of collocation points is M = 27. For the next level, the number 
of collocation points is doubled (the value J is increased by one). The new 


values for u(), u'(x),..., u(x) and a, are estimated as 
2M 1 ; 
AMD a) = Soa fener ig Oak, (16) 
p=1 : 
and form =0,...,n—1 
2M 1 ; 1 
AD (a) = Saf Ggymelthg( Gag — a" fie) hy Ed 
oz 0 0 
n—m—-1 1 
+ pe ae (0), (17) 


where 


a(J—-1 
pale ) forp=1,...,M, (1s) 


0 forp= M +1,...,2M. 


These estimates are corrected by Newton method; the application of which 
for! =1,...,2M leads to the equation 
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shy = — K(x, 0 (x), a (x), ..., a (ay) + F(a), (19) 


OK _ OK du OK Ow | OK dul”) 
Oa, Ou Oap Ou! Oay du”) dap 


By solving (19), Aa,, p = 1,...,2M, are calculated and equations ((4), (7), 
and (18) are corrected as follows: 


1 
uD (2) = S> a) | (i€)"e'**h, (E)dE, (20) 


0 


and 


M 1 1 
uD (x) = Sal’ ie (i€) el, (de — 2} | (i) hp(€) df] 
p=1 


nm—mMm— 


>a <a oylet™(Q) = m=0,...,n—1. 


We start the procedure by taking an initial solution of a = =1 and ao) = = 0. 
These estimates are corrected by solving (1) for J = 0. 


Example 3. Consider the the differential equation 
u(x) + 2u(w)u'(a)? = 2°(1+ V20?), (21) 


which has the exact solution we, (x) = x?. 
The obtained results for J = 2 are shown in Table 3 and Figures 3 and 4. 


Example 4. Consider the the differential equation 
xcos xu’ (x) + (2+ 3)u(x) + xu(x)*? = x(1 +sin x) + 3sinz, u(0) = 0, (22) 


which has the exact solution ue,(x) = sinz. The obtained results for J = 2 
are shown in Table 4 and Figures 5 and 6. 
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Figure 3: The comparison of the numerical solution and exact solution for example 3. 
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Figure 4: The absolute errors for example 3. 
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Figure 5: The comparison of the numerical solution and exact solution for Example 4. 


Figure 6: The absolute errors for example 4. 
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Table 3: Exact solution and The new method solution for equation (2D) 


(a /16) Exact solution The new method solution 
1 0.00390 0.00390 
3 0.03515 0.03515 
5 0.09765 0.09765 
7 0.19140 0.19140 
9 0.31640 0.31640 
11 0.47265 0.47266 
13 0.66015 0.66016 
15 0.87890 0.87890 


5 Error analysis 


In this section, the error analysis for the new method has been offered. 


Lemma Let f € L?(R) be a continuous function defined in (0,1). Then 
the error norm at Jth level satisfies the following inequalities 


2 
E;< 3K 4-3 ; 
16 


where f(x) < K, for all x € (0,1) and K > 0, and M is a positive number 
related to the Jth level resolution of the wavelet given by M = 2/ 


Proof. 


Eyl = lular) —uz(a)] =| S> ap | chy (Cdk, 


p=2M+1 


where 
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Table 4: Exact solution, The new method solution for equation (2) 


(x/16) Exact solution The new method solution 
1 0.0624593 0.0624592 
3 0.1864032 0.1864033 
5 0.3074385 0.3074385 
7 0.4236762 0.4236762 
9 0.5333026 0.5333026 
11 0.6346070 0.6346070 
13 0.7260086 0.7260086 
15 0.8060811 0.8060811 


zz |? = | (a | ede) S> ap | ete hy (€)dé)de 


p=2M+1 p=2M+1 
oo oo 1 1 1 
= YY aya [Cf e8hp(eaey f esia(€)ae)ae 
p=2M+1k=2M+41 o 0 0 
1 2M Tv. 2M 


~ ‘ . apa f Qu holes) f dN he(es) | “ e't§dé)dax. 


7 
p=2M+1 k=2M+4+1 se 
Applying mean value theorem 


SY ante [AaB (wae 


p=2M+1k=2M-+1 


where 
aM 4 
A,(x) = S- ag tv (#s) (cos(ans) + isin(xras)) 
g=1 
and 


2M 
1 


By (x) = S- aag te (ee) (cos(ane) — isin(xag)) 
k=1 
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where 7%, @% © [Tk, 741]. Therefore 
Ze ApGk 
I< BEL YS aglesyhalosy 


p=2M+1k=2M+1 s=1 s=1 


Co 


|a,|? 


Vi 


p=2M+1 


Ay = | 29/? hy (€) f(E)de = | 29/? hy (€)Re(f(€))dé +4 | 24/?h,(€)Im(f(€))dé. 
0 0 0 


By applying the mean value theorem for 6; € (£1, 2) and Be € (€2, €3) 


1 ‘: ; &2 7 &3 “ 
| 2/hn(g)Re(fle)yag = 21] Re(Fle))ag— J Re(Fle))as 
= 21/71(& — &)Re(f(1)) — (6 — &)Re(f(42))]- 


Applying the mean value theorem 


i) 9i/2h, (€) Re( F(€))dé = 2-4/2-1(8, — Ba)(Ref (au), 


0 


where A; € (81, 82) and also 
ak 
‘ 25/2hg(€)Lm(F(E))aé = 2-4/2-1(By — Ba) (If (2) 


where A2 Sc (83, Ba), G3 € (€1, €2), and fon € (&2, &3). 


Therefore 
1 
Ap = | 29/* hy (€) f(E)dé = 274/21 (81 — Bo) Ref (A1) + i(83 — Ba) Imf(A2))’. 
This implies that 
leg cio 


Therefore, 


Co Co 


ere yo Se ye 


i=2M+41 i=2M+1 


A New Hybrid Method Based on Pseudo Differential. . . 79 


OS -/. Ke. = — 
eT, yo. a Sag So 2-8 (ait! 1-2) +1) 
g=IH1 i=2941 j=J+1 
2 x 2 9—-2(J+1) 
gi Spee ee 
M jaJ+l M 1-— q 
3K? 535 
= 06 


It is obvious that the error bound is inversely proportional to the level 
of resolution J of Haar wavelet. Hence, the accuracy in our new method 
improves as we increase the level of resolution J. 

Conclusions In this paper a new method for numerical solution of ordi- 
nary differential equations based on pseudo differential operators and Haar 
wavelets is proposed. Its applicability and efficiency are checked on some 
problems. It is shown that with a small number of collocation points, more 
accurate solutions in comparison with Haar wavelet method can be obtained. 
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